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Abstract 

The ATLAS and CMS experiments did not find evidence for Super- 
symmetry using close to 5/fb of published LHC data at a centre- of -mass 
energy of 1 Te V. We combine these LHC data with data on B® — > n + n~ 
(LHCb experiment), the relic density (WMAP and other cosmological 
data) and upper limits on the dark matter scattering cross sections on 
nuclei (XENON100 data). The excluded regions in the constrained Min- 
imal Supersymmetric SM ( CMSSM) lead to gluinos excluded below 1210 
GeV and dark matter candidates below 220 GeV for values of the scalar 
masses (mo) below 1500 GeV. For large mo values the limits of the 
gluinos and the dark matter candidate are reduced to 910 GeV and 130 
GeV, respectively. If a Higgs mass of 125 GeV is imposed on the fit, the 
preferred SUSY region is above this excluded region, but the size of the 
preferred region is strongly dependent on the assumed theoretical error. 

1 Introduction 

Supersymmetry (SUSY) is a good candidate for physics beyond the SM, because it solves the hierarchy 
problem, allows for unification of the coupling constants and predicts electroweak symmetry breaking 
(EWSB) with the lightest Higgs mass below 130 GeV (for reviews, see e.g. pH]). In addition the 
Lightest Supersymmetric Particle (LSP) has all the properties expected for the Weakly Interacting 
Massive Particles (WIMPS) of the dark matter [SH7], which is known to make up more than 80% of 
the matter in the universe [8]. From the observed DM density, which is proportional to the inverse 
of the annihilation cross section one finds a strong constraint, if one assumes all DM is made from 
supersymmetric LSPs. Unfortunately, direct searches for the predicted SUSY particles at the LHC 
running at 7 TeV were unsuccessful. Also direct DM searches in deep underground experiments were 
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Table 1: List of all constraints used in the fit to determine the excluded region of the CMSSM 
parameter space. 

contradictory [7J. Combining all data from the LHC, cosmology and direct DM searches leads to 
strong constraints on the masses of the predicted SUSY masses, as discussed in many recent papers 

PHZD]- 

To restrict the number of independent SUSY masses one usually assumes the masses to be unified at 
the GUT scale and the particles get different masses at lower energies because of radiative corrections. 
This works well for the SM model particles of the third generation, if they are in the same multiplet 
and get mass from the same Higgs field. E.g. the ratio of b/r masses is well predicted by radiative 
corrections, if one assumes the Yukawa couplings are unified at the GUT scale. For the mass breaking 
terms of the SUSY particles one assumes that the masses of spin (spin 1/2) particles are unified 
at the GUT scale with values w-o( m -i/2)- I n the constrained Minimal Supersymmetric SM (CMSSM) 
|2H [2"2~] the many parameters of SUSY models are reduced to only four: the two mass parameters 
mo, mi/ 2 and two parameters related to the Higgs sector: the trilinear coupling at the GUT scale 
Aq, and tan/3, the ratio of the vacuum expectation values of the two neutral components of the two 
Higgs doublets. Electroweak symmetry breaking (E WSB) fixes the scale of /i, so only its sign is a free 
parameter. The positive sign is taken, as suggested by the small deviation of the SM prediction from 
the muon anomalous moment, see e.g. |23j . 

In this letter we combine the data from LHC, WMAP, XENON100, flavour physics and g-2. 
The specific observables are detailed in Table [T] We start by discussing the fitting technique, the 
observations and the excluded regions of each observation separately. 

2 Multistep Fitting Technique 

Excluded regions have been determined by many different groups either using a frequentist approach 
by maximizing a likelihood or using random sampling techniques of the parameter space, see e.g. 
[9TI20] and references therein. Bayesian techniques, as typically used with Markov Chain or Multinest 
sampling techniques, are dependent on the prior, which leads to an additional, non-quantifiable uncer- 
tainty in the excluded or allowed regions, see e.g. [21] for a recent discussion and references therein. 
We believe this uncertainty is due to the high correlations between three of the four parameters, as 
we discussed in two previous papers [23 . Such strong correlations lead to likelihood spikes in the 
parameter region, where three of the four parameters have to have specific correlated values. Although 
the likelihood of such narrow regions is high, they can be easily missed in methods based on stepping 
techniques. 

To cope with the strong correlations we use a multistep fitting technique, defined by fitting the 
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parameters with the strongest correlation first, i.e. we fit first tan/3 and Aq for each pair of the mass 
parameters mo and rai^ by minimizing the x 2 with the program Minuit [36]. The most probable 
region of parameter space is determined by the minimum Xmin vame - F° r t wo degrees of freedom the 
given 95% C.L. (90% C.L.) limit is reached for an increase in x 2 of 5.99 (4.61). The \ 2 function is 
defined as: 

(nh 2 -0.1131) 2 (&^X,7-3.55-l(T 4 ) 2 



X 2 '2 

(B u -try- 1.68 • 10" 4 ) 2 (Aa M - 302 • 10~ n ) 



+ 



+XsO^. At + M - + Xm h + XCMS + XATLAS + Xm A + XDDMS 

For the first 4 terms in eq. [l]the x 2 is defined in a straight forward way: the square of the difference 
between the predicted value for given SUSY parameters and the experimental value, weighted by the 
inverse of the error squared. For the remaining terms only 95% C.L. or 90% C.L. limits have been 
given. All experimental values and limits have been summarized in Table [TJ We discuss the remaining 
contributions to the x 2 separately: 

• X 2 B o^+^- = {B° s fM+fJ,- - B° s -> fi + fi 95 ) 2 /a% 

It is added only if B® — > fi + pT > B® — > [i + fig 5 . B® — > /^ + ^g 5 was set such that with the 
theoretical error ob = 0.2 • 10~ 9 [29J a change in x 2 of 5.99 was reached for the experimental 
95% C.L. upper limit from Table [lj This leads to B° s -»■ = 4.0 • 10" 9 . 

• Xm h = {m h -mf) 2 /al 

This contribution takes the lower limit on the Higgs mass into account. It is added to the x 2 
function if m,h < rrv^. Here a 95% C.L. lower limit of rrih of 114.4 GeV is used, as determined 
from the LEP experiments. The experimental error taken from the la band from Ref. |30| 
corresponds to the chosen Oh of 0.5 GeV, which result in m 95 = 115.6 GeV from the requirement 
A X 2 =5.99 for m h > 114.4 GeV at 95% C.L.. 

• Xcms an d Xatlas 

These contributions take into account the negative results from the direct searches for SUSY 
particles from CMS and ATLAS, respectively. Low SUSY masses are excluded by these searches, 
since these would have a too large cross section. From the published 95% C.L. exclusion contours 
in the (mo, m 1 / 2 )-plane [33j[3l] one can determine the excluded SUSY cross sections, which vary 
along the contour presumably because of the varying efficiency. Since the efficiency is not 
published, we take the following simplified approach for the x 2 contribution: we know that at 
the contour A% 2 =5.99 and we assume that the exclusion limit near the contour is proportional to 
the SUSY hadronic cross section a tot (pp — > gg, gq, qq), which is a reasonable approximation for 



the LHC limits as will be shown later in section 3.2 Hence we parametrize Xlhc = a tot/ a LHC> 
where LHC stands either for CMS or ATLAS and for each experiment <?\ HC was determined as 
function of tuq and m^2 by the requirement that A% 2 =5.99 at the 95% C.L. contour. 

v 2 

This contribution takes care of the results of the search for the neutral Higgs bosons decaying 
into tau pairs in pp collisions at the LHC. The 95% C.L. exclusion curve in the tan/3 — plane 
excludes for a given value of tan (3 a certain mass of the pseudo-scalar Higgs [311 [32] . The x 2 
contour is calculated as x 2 = (m/^ (tan /3) — m 9 4 5 (tan/3)) 2 /o"^g 5 , where c m 95 can be obtained from 

the la curves around the 95% C.L. contour line. m 9 4 5 (tan/3) is determined by the requirement 
that Xrn A 1S 5.99, for m^(tan/3) on the contour line. 
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Figure 1: Left: A% 2 distribution of all constraints up to m^ = 3000 GeV. The white cross represents 
Xmin = 4.06. The white region in the top left corner is excluded because the stau is the 
LSP. The red region in this corner is excluded by the relic density constraint requiring large 
tan (3, which in turn causes a large mixing in the stau sector leading to the stau becoming 
the LSP again. Right: contributions to the x 2 of all contraints up to mi/ 2 = 1000 GeV. 
The contour for each constraint represents the 95% C.L. exclusion limit (Ax 2 = 5.99) for 
each constraint separately. 



This term takes the upper limits on the elastic WIMP-nucleon cross section <j x n of the direct 
dark matter searches into account. Usually only 90% C.L. upper limit on a x jy are given as 
function of the WIMP mass mwiM p [35] . The x 2 contribution can be included in a way similar 
to Xm A - 

We define x x n = ( CT YJv( m x) ~~ ^viv^x)) 2 V ,<7 xiv( m x) ■ ^he we ight l/c 2 ^ is taken from the la 
band in Ref. [35J. The excluded cross section cr x 5 N (m x ) is determined by requiring \ x n =4-61, 
since the published limit is at 90% C.L.. 

All observables were calculated with the public code micrOMEGAs 2.4.5 [371 EE] combined with 
Suspect 2.41 as mass spectrum calculator [39]. Within our multistep fitting technique we minimizie 
the x 2 function defined in eq. [TJ which is fast, since initially only two parameters are fitted for each 
point in the (mo,Jn 1 / 2 )-grid. 

Once the SUSY parameters have been fitted, one should also vary the SM parameters or marginalize 
over them, like the top and bottom mass and the strong coupling constant. The SM parameters are 
given in the Particle Data Book [40J: we use m v °^ = 172.5 ± 1.3 and m b (m b ) MS = 4.25 ± 0.2 GeV 
for the heavy quark masses and a s = 0.1172 ± 0.02 for the strong coupling constant. However, the 
SM quark masses mainly determine the running of the Higgs mass parameters, so different values of 
these masses can be compensated by a slightly different values of the SUSY parameters Aq and tan f3 
in order to get the same x 2 value, so the allowed region is hardly affected. 

3 Combination of all Constraints 

The minimal values of the x 2 function of Eq. [I] in the (mo, m 1/ / 2 )-plane are shown in the left panel 
of Fig. [ljfor each mo,m]/ 2 pair. The colour coding indicates the value of A^ 2 = X 2 ~ Xmin' where 
the the minimal value Xmin = 4-06 is obtained for the mass values mo = 350 GeV, m 1( / 2 = 825 
GeV, as indicated by the white cross. The red and yellow regions, corresponding to A^ 2 > 5.99, are 
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Figure 2: Corresponding tan/3 (left) and Aq (right) values to the fit to all data. 

excluded at 95% C.L.. The observables contributing most strongly to the exclusion vary in the plane 
as indicated in the right panel of Fig. [TJ These contours are drawn in the following way: we take 
A-Xi = xj ~ Ximin = 5.99, where Ximin ^ s * ne X 2 - con t r ibution of variable i at the best-fit point and 
xj is the x 2 value of variable i at the contour. The direct SUSY searches at the LHC (contour 1) 
dominate the limit at small tuq with a small contribution from the branching ratio of B® — > n + \i~ 
(contour 2). Other contributions at intermediate SUSY masses come from the Higgs searches (contour 
3 for the SM Higgs and contour 4 for the pseudo-scalar Higgs) and direct DM searches (contour 5) at 
larger values of tuq. The fitted values of tan /3 and Aq for each mo, mi/2 pair are shown in Fig. [2j 

To understand the contours in the right panel of Fig. [TJ we discuss each of them in more detail 
after discussing the influence of g-2 first. 

3.1 Influence of the anomalous magnetic moment of the muon 

We included the value of g-2 into the fit, but as can be seen from Fig. [3]most of the preferred region by 
the g-2 constraint is excluded by the direct searches of the LHC independent of the treatment of the 
g-2 uncertainties. The green region (dark shaded) is preferred by g-2 data if one adds the experimental 
and theoretical errors in quadrature. However, since these errors are of the same order of magnitude 
(see Table. [TJ and the theoretical uncertainties are certainly non-gaussian, a linear addition of the 
errors is more conservative, which leads to a larger error and hence a larger (yellow, light shaded) 
preferred region. However, even this larger region is still excluded by the LHC, so the observed three 
sigma deviation of the anomalous magnetic moment of the muon above the SM prediction may either 
be a statistical fluctuation or has an origin different from SUSY, if we assume the theoretical and 
experimental errors have been estimated correctly. 

3.2 Excluded region by direct searches for SUSY at the LHC 

The direct searches for Supersymmetry at the LHC are dominated by the search for strongly interacting 
particles, as shown in the publications by ATLAS |33j and CMS [HI]. From Fig.[4]one observes that the 
excluded region (below the solid line) follows rather closely the total cross section for the production 
of squarks and gluinos atot, indicated by the colour shading. The 95% C.L. on a tot {pp — > 9£l>Mi<j<j) 
is given by the contour line and varies between 0.003 and 0.03 pb, as shown in Table [lj This is the 
motivation to approximate the x 2 contribution from the LHC experiments near the contour as a^ ot /a 2 , 
where a is defined by the requirement that at the exclusion limit the x 2 contribution equals 5.99 for 
each LHC experiment, as discussed in section [2] 
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Figure 3: Preferred region of the g-2 observable alone under the constraint that tan/3 and Aq are 
fixed by all other constraints. Here we show the la band (A% 2 = 2.3) of the preferred 
region for quadratic (green, dark shaded) and linear (yellow, light shaded) addition of the 
theoretical and experimental errors. We compare these bands with the 68% C.L. exclusion 
limit of the direct searches at the LHC. The preferred region by g-2 is largely excluded by 
the LHC constraints. 



If the LHC data is combined with cosmological and electroweak data, the fitted values of tan /3 and 
the trilinear coupling vary, while the LHC experiments provided the exclusion contours for fixed values 
tan/3 = 10 and Aq = 0. However, the efficiencies are not very sensitive to these parameters, so 
this definition of Xlhc nas ^ ne advantage that the exclusion curve for each experiment is perfectly 
reproduced and that the dependence of the cross section on Aq and tan/3 is taken into account. By 
adding the exclusion from each LHC experiment independently to the x 2 we assume tacitly that 
there are no correlations between these independent data sets. If we assume a 100% correlation, we 
should only use a single experiment. However, the results would hardly change, as can be seen from a 
comparison from the combined curve in Fig. [I] and the published curves from each of the experiments: 
the difference in m^ 2 is at most 25 GeV. 

3.3 Excluded region by B® — > 

The upper limit on the branching ratio of B® — > /i + can give significant constraints on the SUSY 
parameter space, since the B® — > rate varies as tan 6 (3. In addition B® — > is sensitive to 

the stop mixing which is a function of Aq. B® — > /z + [i~ can be suppressed if m t ~ ~ mf or even get 
values below SM |25j . Hence the \ 2 is sensitive to the chosen tan (3 and Aq value. The combination 
with the relic density, which requires a large tan (5 value in a large region of parameter space (see Fig. 
[2J causes tension with the B® — > /x + /i~ constraint. This tension can be reduced by large values of Aq 
as we showed in |25| . but with the recent upper limit near the SM value from LHCb [29J this tension 
increased and both constraints cannot be fulfilled at the same time in the whole parameter space. 
This leads to two excluded regions shown in Fig. [l] by contour 2. The reason for the two regions is the 
following: at small values of uiq the trilinear coupling Aq cannot be made large enough to suppress 
B® — > jjL + fj,~ enough, because then the stau leptons become tachyonic. At intermediate values of too 
the trilinear couplings can be made large enough to suppress B® — > /i + /U _ , but at larger values of tuq 
rriA becomes too large for large Aq values, which leads to a too large relic density. Compared to the 
other constraints the excluded region by B® — > /i + /i _ leads only to a tiny increase of the excluded 
region at small rriQ. 
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Figure 4: Left: Total production cross section of strongly interacting particles (colour coding) in 
comparison with the LHC excluded limits for 7 TeV. Here the data from ATLAS and 
CMS were combined. The ATLAS and CMS data correspond to an integrated luminosity 
of 4.4 and 4.71 fb _1 , respectively. One observes from the colour coding that a cross section 
of 0.003 to 0.03 pb is excluded at 95% confidence level. Right: Values of tua in the 
(mo, mi/2)-pla n e after optimizing tan/3 and Aq to fulfill all constraints at every point. 
The data below the solid line in the right panel is excluded at 95% confidence level from 
the rriA exclusion contour as function of tan /3. 




Figure 5: If a Higgs mass of 125 GeV is imposed in the fit, the best-fit point moves to higher SUSY 
masses, but the location is strongly dependent on the assumed error for the calculated 
Higgs mass. This error is indicated by the number inside the circle for the best-fit point. 
Left Ax 2 = 2.3(lcr) contour; right Ax 2 = 5.99(2cr) contour. 



3.4 Effect of a SM Higgs mass rrih around 125 GeV 

The 95% C.L. LEP limit of 114.4 GeV contributes for small and intermediate SUSY masses to the % 2 
function, as shown by contour 3 in Fig. [U In the fit we use the 95% C.L. LEP limit of 114.4 GeV 
on the Higgs mass instead of the limits published by CMS and ATLAS with about 5/fb. In these 
publications CMS [41] and ATLAS |42] show some evidence for a Higgs with a mass around 125 GeV. 
If we assume this to be evidence for a SM Higgs boson, which has similar properties as the lightest 
SUSY Higgs boson in the decoupling regime, we can check the consequences in the CMSSM we are 
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investigating. If a Higgs mass of 125 GeV is included in the fit, the best fit point moves to higher 
SUSY masses, but there is a rather strong tension between the relic density constraint, B® — > fi + [i~ 
and the Higgs mass, so the best fit point depends strongly on the error assigned to the Higgs mass, 
as shown in Fig. [5] The experimental error on the Higgs mass is about 2 GeV, but the theoretical 
error can be easily 3 GeV. Therefore we have plotted the best fit point for Higgs uncertainties between 
2 and 5 GeV. One sees that the best fit points wanders by several TeV. Clearly this needs a more 
detailed investigation in the future. It should be noted that the fit does not provide the maximum 
mixing scenario. If we exclude all other constraints, the maximum value of the Higgs can reach 125 
GeV, albeit also at similarly large values of m^ 2 . A negative sign of the mixing parameter [i shows 
similar results. 



3.5 Excluded region by the pseudo-scalar Higgs mass tua 

The pseudo-scalar Higgs boson mass is determined by the relic density constraint, because the dom- 
inant neutralino annihilation contribution comes from A-boson exchange in the region outside the 
small co-annihilation regions. One expects itla oc m^2 from the relic density constraint, which can 
be fulfilled with tan/3 values around 50 in the whole (mo, mi/2)-plane [26]. Since the vtia production 
cross section at the LHC is proportional to tan 2 (3 the pseudo-scalar mass limit increases up to 496 
GeV for the large values of tan/3 preferred by the relic density (see Fig. [2]). 

In our fit we are not only using the relic density as a constraint but all data. The result of this 
optimization was shown in Fig. pleading to different values of Aq and tan/3 in the (mo, m 1( / 2 )-plane. 
The corresponding m^-values are displayed in the right panel of Fig. [4] and the rriA values excluded 
by the LHC searches lead to the excluded region, shown by the contour line in Fig. [4] (identical to 
contour 4 in Fig. [T] right). 

3.6 Excluded region by direct DM searches 

The cross section for direct scattering of WIMPS on nuclei has an experimental upper limit of about 
10 -8 pb, i.e. many orders of magnitude below the annihilation cross section. This cross section is 
related to the annihilation cross section by similar Feynman diagrams. The many orders of magnitude 
are naturally explained in Supersymmetry by the fact that both cross sections are dominated by Higgs 
exchange and the fact that the Yukawa couplings to the valence quarks in the proton or neutron are 
negligible. Most of the scattering cross section comes from the heavier sea-quarks. However, the 
density of these virtual quarks inside the nuclei is small, which is one of the reasons for the small 
elastic scattering cross section. In addition, the momentum transfer in elastic scattering is small, so 
the propagator leads to a cross section inversely proportional to the fourth power of the Higgs mass. 

Since the particle which mediates the scattering is typically much heavier than the momentum 
transfer, the scattering can be written in terms of an effective coupling, which can be determined 
either phenomenologically from irN scattering or from lattice QCD calculations. The default values 
of the effective couplings in micrOMEGAs [43] are: = 0.033, /| p j = 0.023, = 0.26, f$ = 
0.042, fftj = 0.018, f$ = 0.26. The lower values from the lattice calculations Q] are: = 

0.020, fil = 0.026, = 0.02, f$ = 0.014, f!$ = 0.036, f!$ = 0.02. Hence the most important 
coupling to the strange quarks vary from 0.26 to 0.02, which implies an order of magnitude uncertainty 
in the elastic neutralino-nucleon scattering cross section. 

Another normalization uncertainty in direct dark matter experiments arises from the uncertainty 
in the local DM density, which can take values between 0.3 and 1.3 GeV/cm 3 , as determined from the 
rotation curve of the Milky Way, see Ref. |45} 06] and references therein. 

To get conservative estimates for the excluded regions, we take the lowest possible values of the 
local DM density and the low couplings from lattice QCD calculations. The excluded region from the 
XENON100 cross section limit |35| is shown by the contour line 5 in Fig. [TJ At large values of mo 
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EWSB forces the higgsino component of the WIMP to increase and consequently the exchange via 
the Higgs, which has an amplitude proportional to the bino-higgsino mixing, starts to increase. This 
leads to an increase in the excluded region at large mo and has here a similar sensitivity as the LHC. 
If we would take the less conservative effective couplings from the default values of micrOMEGAs the 
XENON100 limit would be 50% higher than the LHC limit. 

4 Summary 

As mentioned in the introduction, several groups have performed similar analysis. Our results are 
closest to the one of Ref. |10j if we adjust to the lower luminosity of 1/fb luminosity used in their 
analysis. However, in Ref. [10] values of m\/2 above 2500 GeV are excluded due to the tension with 
the relic density constraint [U]. In our case we do find good solutions and no such excluded region is 
found as shown in Fig. [TJ left panel. This is probably due to the fact that in this region tan f3 and Aq 
are highly correlated, so they can be easily missed, if SUSY samples are prepared with many, but not 
necessarily all, combinations of parameters. This region is the one preferred for Higgs masses around 
125 GeV, as shown in Fig. [5j 

Our results differ significantly from results using Markov Chain Monte Carlo sampling. E.g. in 
Ref. [11] values for intermediate values of mo are excluded, which is the region of large tan/3 (see 
Fig. [2j left panel). Here the parameters tan/3 and Aq are highly correlated and finding the correct 
minimum depends strongly on the stepping algorithm, e.g. stepping in the logarithm of a parameter 
is different from stepping in the parameter ("prior dependence"), see e.g. Ref. [17] . Such dependence 
on sampling techniques largely disappears in our multistep fitting technique, since for each point of 
the (mo, mi/2)-g r id a unique solution is found independent of the minimizer used, so the frequentist 
approach with x 2 minimization yields the same results as a likelihood optimization with a Markov 
Chain sampling technique. 

If one combines the limits from the direct searches at the LHC, heavy flavour constraints, WMAP 
and XENON100 we exclude values of m 1/2 below 525 GeV in the CMSSM for m < 1500 GeV, which 
implies a lower limit on the WIMP mass of 220 GeV and a gluino mass of 1270 GeV, respectively. For 
larger values of mo the excluded region drops to mi/2 below 350 GeV, which leads to a lower limit on 
the LSP mass of 130 GeV and a gluino mass of 970 GeV, respectively. 

If a Higgs mass of the lightest Higgs boson of 125 GeV is imposed, the preferred region is well 
above this excluded region, but the size of the preferred region is strongly dependent on the size of the 
assumed theoretical uncertainty, as shown in Fig. [5] However, the higher than expected branching 
ratio and lower than expected branching ratio into tau leptons point to a Higgs with slightly different 
couplings from the SM. Such different couplings could exist in a supersymmetric model with an 
extended Higgs sector, like the NMSSM, see e.g. @8J. 
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6 Note added in proof 

The CMS and ATLAS collaborations confirmed the published evidence for a new boson with a mass 
around 126 GeV at the beginning of the ICHEP 2012 conference [3§], so the CMSSM would indeed 
prefer heavy SUSY masses, as shown in Fig. [5j This conclusion does not change, if one relaxes the 
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CMSSM constraints, e.g. by not requiring unification of the Higgs masses at the GUT scale. However, 
extended Higgs sectors would need additional investigation. 
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